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Although lower hybrid waves are effective at driving currents in present-day tokamaks, they are expected to 
interact strongly with high-energy particles in extrapolating to reactors. In the presence of a radial alpha 
particle birth gradient, this interaction can take the form of wave amplification rather than damping. While 
it is known that this amplification more easily occurs when launching from the tokamak high-field side, the 
extent of this amplification has not been made quantitative. Here, by tracing rays launched from the high- 
field-side of a tokamak, the required radial gradients to achieve amplification are calculated for a temperature 
and density regime consistent with a hot-ion-mode fusion reactor. These simulations, while valid only in 
the linear regime of wave amplification, nonetheless illustrate the possibilities for wave amplification using 
high-field launch of the lower hybrid wave. 

PACS numbers: 52.35.-g, 52.55.Fa, 52.55.Wq, 52.55.-s 


I. INTRODUCTION 

The lower hybrid (LH) wave has been show n to be 
effective in driving plasma currents in tokamaks.l^ Yet 
there remains a concern that in a fusion reactor, high- 
energy a particles born in the plasma core could strongly 
damp the LH wave.^^^ However, it is possible for a favor¬ 
able population inversion to appear along the diffusion 
path, resulting in wave amplification rather than wave 
damping.^ This amplification effect, which relies on cou¬ 
pling diffusion in energy to diffusion in space, is often 
referred to as alpha channeling. 

For lower hybrid waves, this amplification appears to 
be best achieved by launching the wave from the tokamak 
high-field side.^ Such a launch is often referred to as “in¬ 
side launch,” since the magnetic field is strongest near the 
tokamak hole. Inside launch represents an engineering 
challenge, since there is less space near the tokamak hole 
in which to place waveguides. However, there are also 
several recently-discovered advantages to inside launch, 
even apart from the opportunity for a channeling, includ¬ 
ing deeper plasma penetration of the LH wave and a m ore 
protected plasma environment for the LH waveguides .1^ 
When it is launched from the high-field side, the oppor¬ 
tunity to amplify the LH wave through interactions with 
a particles adds to these advantages. However, while this 
amplification indeed appears to more easily occur when 
launching from the tokamak high-field side, the extent of 
this amplification has not yet been made quantitative. 

To quantify the conditions under which ray amplifi¬ 
cation can occur, it necessary to consider the ray paths 
for the lower hybrid waves, including both damping on 
electrons and growth on a particles along the ray tra¬ 
jectory. Of course, nonlinear effects can affect impor¬ 
tantly the propagation and damping of the lower hybrid 
waves, particularly with respect to penetration of the hot 
and dense core of the tokamak reactor. These nonlinear 
effects include quasilinear flattening of the distribution 
function which would lessen the electron dampin^iSl and 


parametric effects that might allow the wave to access 
higher density.^ The employment of the linear model 
here is thus not meant to optimize penetration to the 
plasma core; rather, it is meant to clarify the new phys¬ 
ical effects in the simplest model that captures these ef¬ 
fects. It should be noted that the amplification effects 
due to alpha channeling, to the extent that linear damp¬ 
ing by electrons is opposed, will in fact enable deeper 
penetration of the wave. In a given tokamak magnetic 
geometry, with given density and temperature profiles, 
the ray paths may be calculated using GENRAY.E^ For 
simplicity, near-circular, near-concentric flux surfaces are 
assumed. Once the ray trajectories are known, both the 
damping and the growth can be separately calculated. 

In the linear approximation for the waves, which for 
simplicity we adopt here, the effect of the LH waves on 
the distribution functions of electrons, ions, or a parti¬ 
cles is ignored. Generally in the regimes of interest here, 
damping on fuel ions is neglected. In terms of the elec¬ 
tron interaction, our linear approximation means that we 
ignore the formation of a quasilinear plateau, which could 
lessen the electron damping. We similarly ignore the re¬ 
laxation of the a-particle distribution, which would limit 
the energy extractable from the a particles. This linear 
approximation is thus valid in the limit of low power in 
the LH wave, where diffusion by waves is small compared 
to diffusion by collisions, so that the ions and electrons 
assume Maxwellian velocity distributions, while the high- 
energy a particles assume a slowing-down distribution. 

Thus, although limited to low LH wave power, the lin¬ 
ear model is used to demonstrate how inside launch can 
lead to dramatic a-mediated LH wave amplification in 
the reactor regime at reasonable radial a-particle birth 
gradients. The very important generalization of these 
effects to higher powers, or for that matter to more re¬ 
alistic tokamak magnetic configurations, is outside the 
scope of this work. Instead, we demonstrate that the lin¬ 
ear damping rate on electrons may be exceeded by the 
linear growth rate on a particles, so that the wave under¬ 
goes significant convective amplification before it finally 
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is extinguished by linear damping. Apart from its own 
intrinsic interest, the calculation of amplihcation effects 
in the low-power limit gives important insights to con¬ 
ditions for amplihcation in the regime of high-power LH 
waves, which is the regime where the full potential of the 
amplihcation effect would be realized. 

The paper is organized as follows: In Sec. [Ill th® linear 
model is explained in greater detail. In Sec. |HI[ an exam¬ 
ple of joint optimization of current drive and a channeling 
is given. In Sec. |IV[ conclusions and caveats with respect 
to regimes of applicability of the results are discussed. 


II. LINEAR MODEL 


Suppose, for simplicity, a tokamak, with circular, con¬ 
centric hux surfaces. Such a magnetic conhguration is 
conveniently parameterized by the toroidal angle 4>, the 
poloidal angle 0, and the minor radius r; then the hux 
surface normal vector coincides with the minor radius 
vector f. Since a particle density and temperature tend 
to be constant over each hux surface, a particle gradients 
will thus also point along the minor radius vector. 

Alpha channeling occurs when a wave with frequency w 
and wavevector k encounters an a particle orbiting with 
velocity u_l around a magnetic held B. This interaction 
couples particle diffusion in energy and space, such that 
particles that gain energy are pushed in the direction of 
X = k X B, while those that lose energy are pushed to¬ 
wards —X. In a reactor, it is desirable to have x point 
inward along the hux surface normal vector, so that par¬ 
ticles lose energy as they are diffused towards the plasma 
periphery. Consider now the scalar quantity dehned 
consistently with previous worlPl, which encodes infor¬ 
mation about the direction of channeling: 

k X B .. kg 

^“IkxBl'”''^ Ifcxl’ 


where kg is the 6 component of the wavenumber, and 
fc_L is the component perpendicular to the magnetic held. 
When ^ is negative, a particles that lose energy will be 
channeled outwards towards the plasma periphery, as de¬ 
sired. However, if |^| 1, the particles are channeled 

largely poloidally (in the direction of 9 rather than ra¬ 
dially, so that there is unlikely to be a steep a particle 
gradient along the direction of channeling. Thus ideally, 

e«-l. 

This coupling of energetic and spatial diffusion results 
in diffusion paths in r — e_L space, where e_L is the a par¬ 
ticle kinetic energy perpendicular to the magnetic held.l^ 
For simplicity, we normalize e_L to the a-particle birth 
energy « 3.5 MeV, so that e± = {v±/va)‘^, where Va 
is the a particle birth energy. These diffusion paths obey: 

de± rUai^aUi 2ZaeBuj 

dx kj_ea k±maV^ ’ 


where = 2 is the ion number, e is the electron charge, 
and rria is the a particle mass. Since the density may be 


assumed to be uniformly distributed along the flux sur¬ 
face, only distances along the radial coordinate matter, 
so that, using Eq. 0 , we find 

9e_L 1 i9e_L rriaflaUJ 2ZaeBuj 

dr ^ dx ik^ea ik^rriavl' 

Thus, as ^ gets smaller, i.e. as the a particles are pushed 
by the wave more in the direction of 0, it takes a very 
large change in energy to move a particles even slightly 
radially. 

To calculate the change in power of the wave due to 
the energetic diffusion of the a particles, consider an in¬ 
finitesimal length I of the propagating ray, which we take 
to have cross-sectional area A. Particles will be pushed 
in the direction x, which, since it must be normal to k, 
necessarily lies in the cross-sectional plane; we denote 
the remaining orthogonal direction y. Thus, the total 
number of a particles transferred across a wave section 
of length dl in a unit time dt is = J{ydl)dt, where 
J is the energetic flux, given from the energetic diffusion 
coefficient and the a particle distribution f{e±,x) by 


J = £>, 


d 

de 



d 

dx 




(4) 


As these particles diffuse, they gain or lose (unnormal¬ 
ized) energy Aej_ = ea{de±/dx)x. Multiplying the num¬ 
ber of particles transferred by the energy lost per particle, 
and integrating across the energetic distribution, we then 
have (noting A = xy ): 


d^ 

Id 


de± eaA 





1) 


(5) 

To make use of this equation, we must find expressions 
for the distribution function /, diffusion coefficient D^, 
and area A. 

To calculate the distribution function, note that the 
energetic distribution of high-energy a particles on each 
flux surface is dominated by collisional slowing down on 
electrons. Suppose a flux-surface-dependent birth con¬ 
centration from fusion na{r); suppose furthermore that 
this birth distribution slows down on electrons, in such 
a way that that the a particles remain on their birth 
flux surface. Then the slowing down distribution is in 
perpendicular energy velocity space can be put as 


ff \ '^a(r) VI - e-L s 

f{f^±,r) = ^ - ^ (6) 

where ly = 16\/27rmj(e^ne x InA/STe^^nia is the colli¬ 
sional damping on electrons and B(x) is the Heaviside 
function. As discussed above, throughout this paper, the 
LH wave intensities are assumed insufficient to modify 
the a-particle distribution, so that Eq. 0 holds: we will 
later discuss the consequences of this assumption. 

The diffusion coeffi cien t in perpendicular velocity 
space can be written a^^E^ 


De(e±) 


2uj 


4(e-L - £w) 


1/2 


'^OSC 

Va 


2 

ff(e± - Cu,), 


(7) 
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where Vosc = 2eE/maUj. Note that if 

the perpendicular phase velocity Vp^±_ = a;/fc_L of the wave 
is greater than the birth velocity Va of the a particles, 
then Cu, > 1. Since e_L < 1, the diffusion coefficient 
will then be strictly zero, and there will be no a particle 
interaction. 

Finally, the power in the wave at a point along the ray 
trajectory is given approximately by 

P = \vge^E^A, (8) 


where A is the ray’s cross-sectional area, E is the electric 
held, eo is the permittivity constant, and Vg is the group 
velocity. 

Plugging the expressions for the diffusion coefficient 
(Eq. 7), distribution function (Eq. [^, and ray power 
(Eq. 8) into Eq. ([^ allows us to calculate the change 
in power along the ray due to the a particle interaction: 


dP 

H 


= P 


/ Zle^Tig \ \ / 1 - e^, 

\veomujelvg) ^ ^ V 



-I- TT ^1 -I- 




( 9 ) 


The two terms within the square brackets have intuitive 
physical interpretations. The hrst term, which is strictly 
negative, arises from the monotonically decreasing en¬ 
ergetic distribution of a particles on each hux surface, 
which can only result in wave damping. The second term, 
arising from the radial a particle birth gradient, can take 
either sign. In particular, if the second term is positive 
and larger in magnitude than the first term, the wave 
will exponentially amplify on the a particles. 

Let us define the radial a-particle birth decay length 
Ao = fig/{dfia/dr) at which this transition to amplifica¬ 
tion occurs. Substituting Eq. (|^ into Eq. and noting 
df /dr ~ drig/dr, yields 


Ao — 4 (1 — ^/e//) 



$ 


oc (1 


( 10 ) 


When Ao is negative and the radial decay length Aq, in 
the a-particle birth distribution satisfies Aq < Ao, then 
the a particles will amplify, rather than damp, the wave. 


III. EXAMPLE OF JOINT OPTIMIZATION 


= 2.693. The temperature and density profiles for elec¬ 
trons and ions are given analytically on each flux surface 
by 

T{r) = (T^ax - Tmin) (l " + Ifnin (H) 

Tl(r) = (rZniax ^min) (l ) T^min; (1^) 

where Tmax — 10 keV, Tniin — 500 eV, Umax = 5 X 
10^^ cm“^, and nmin = 10^^ cm“^. Here the radial coor¬ 
dinate r refers to the radius of each (essentially) circular 
flux surface with respect to its own center, not with re¬ 
spect to the magnetic axis. Thus, the density and tem¬ 
perature are constant on a flux surface, where the larger 
radii surfaces have the smallest densities and tempera¬ 
tures. The flux surfaces for this equilibrium exhibit just 
a small Shafranov shift (Ai? < 15 cm), and, while not 
concentric, are essentially circular. Deuterium and tri¬ 
tium ions are in equal proportion. The pressure of a 
particles is neglected. 

With this equilibrium and these kinetic profiles, we cal¬ 
culate one-pass ray trajectories in GENRAY,E^ a geomet¬ 
rical optics code that also calculates the linear Landau 
damping on electrons (Landau damp ing on the ion dis¬ 
tribution is assumed to be negligible). b^b6 | Because geo¬ 
metrical optics assumes the WKB (short-wavelength) ap¬ 
proximation, the effects of beam widthfi^ and short-scale 



As an example of joint optimization of the current- 
drive and a-channeling effects in the reactor regime, con¬ 
sider a nearly circular equilibrium,^^ with parameters 
similar to the proposed ARC tokamak reactor.l^ Accord¬ 
ingly, take major radius R — 3.45m and minor radius 
a = 1.13m, with a toroidal magnetic field of 9.25 T 
and 7.8 MA of plasma current. For this equilibrium, 
the safety factor q monotonically increases with distance 
from the magnetic axis, from g(axis) = 1.015 to g(edge) 


FIG. 1. Simulation results for an ARC-like equilibrium for 
ny = —1.65, / = 3.7 GHz, and 6 = -k (inside launch), (a) 
Inset: poloidal trajectory of ray. Graph: the evolution of kg 
(blue), \kr\ (red), and k± (green) along the poloidal trajectory 
length s. Near the reflection point (s ~ 2.1 m), most of fcx is 
in kg, and so ^ becomes large and negative, (b) Evolution of 
As the ray gets near the center, the increase 
in fcx (and thus decrease in ux) causes the wave to interact 
with the hot a particles. Color available online. 
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FIG. 2. The evolution of the components of the group veloc¬ 
ity for an ARC-like equilibrium for = —1.65, / = 3.7 GHz, 
and 6 = -K (inside launch). Gomponents are Vg^g (blue), |ug,r| 

(red), and U 9 ,poi = (gr^en). Near the reflection 

point, the poloidal group velocity is dramatically reduced, so 
that the trajectory spends far “longer” where channeling is 
most favorable. Color available online. 


inhomogeneity on wave propagation are not captured. 
Additionally, the code employs a cold-plasma dispersion 
relation, so potential warm-plasma effecti^^ on the ray 
trajectory are not considered in the current study. How¬ 
ever, both these effects could be incorpor ated b y combin¬ 
ing our linear model with full-wave code^i^f^Sl employing 
a warm-plasma dispersion relation. 

We assume that the ray encounters a constant a- 
particle gradient, arising from an a particle birth distri¬ 
bution with decay length in the region of wave prop¬ 
agation, and say constant otherwise, i.e. 

. J naoexp(-r/Aa) if r > Cdosest 

= oth™ ' 


where r’dosest is the distance of closest approach of the ray 
to the plasma center. The linear damping (or growth) 
rate due to a particles can then be calculated along the 
ray trajectory, taking into account the linked diffusion 
in r — ej^ space, according to Eq. (§. Here we see the 
main simplification of the circular, concentric equilib¬ 
rium, since the local radial coordinate, with Shafranov- 
shifted origin, can then be used as the flux normal vector 
for each surface. 

Consider, for example, a peak high-energy a particle 
density (above e; = 30 keV) of Ua = //“ /(e_L)de_L = 
7.6 X 10^^ cm“^, a radial decay length |Aa| = 0.8 cm, and 
^closest = 8.43 cm. Then consider launching an LH wave 
from 9 = TT (inside launch) with ny = —1.65 and / = 3.7 
GHz. As this ray gets close to its reflection point near the 
plasma center, it experiences a sharp increase in kg and 
a sharp decrease in (Fig. leading to an increase in 
the magnitude of the component of channeling normal 
to the flux surface. Because Ao oc this effect leads 
to far more favorable conditions for a channeling near 
the reflection point. Crucially, this transfer of k into k± 
occurs as kj_ reaches its peak value, causing e^, to drop 
below the interaction threshold at e^, = 1, as seen in 


Fig.[T]3. Thus the wave interacts most strongly with the 
a particles precisely at the point where Aq becomes most 
favorable for channeling. 

Simultaneously, the ray also experiences a dramatic 
decrease in the poloidal group velocity (Fig. [^, causing 
the ray to spend far longer near the plasma center than 
suggested by the distance it travels in the poloidal plane, 
exactly where channeling is most favorable. 

We can see the importance of this combination of fac¬ 
tors in Fig. Near the plasma center, Aq increases dra¬ 
matically, resulting in a transition from wave damping 
to wave amplification as shown in Fig. [3}r. The ray thus 
regains power from the a channeling effect on the order 
of its initial power as shown in Fig. Wp- As a result, 
the absorbed power by the electrons is more than twice 
the launched power in the ray. Furthermore, this excess 
power is all damped on electrons near the plasma center 
(see Fig. |^, which is more favorable for plasma stability. 

It was previously found from essentially analytic con¬ 
siderations that launching from the high-field side, al¬ 
though unconventional, would be critical to achieve the 
wavenumbers needed for a channeling.^ Our present sim¬ 
ulations confirm this finding. As can be seen in Fig. 
wave launch from the low-field side results in large neg¬ 
ative values of kg near the plasma center. This rever¬ 
sal in the poloidal wavenumber means that, although 
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FIG. 3. Simulation results near the plasma center for an 
ARG-like equilibrium for ny = —1.65, / = 3.7 GHz, and 
8 = TT (inside launch), (a) Gyan line: negative radial decay 
length in cm (Ao) of a particles required for zero damping vs. 
poloidal length along ray. As the ray nears the plasma center 
(s « 2.0 m), the maximum allowable decay length increases. 
Black line: imposed radial decay length of Aq = —0.8 cm. (b) 
Power remaining in ray (blue) and cumulative power delivered 
to electrons (red) along ray trajectory, both with (solid) and 
without (dashed) a particles present. Golor available online. 
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FIG. 4. Radial damping profile on electrons for an ARC- 
like equilibrium with ny = —1.65, / = 3.7 GHz, and 0 — tt. 
Dashed line shows the damping profile in the absence of a 
particles, solid in the presence of a particles. Most of the 
power resulting from a. channeling is damped near the plasma 
center. 


Aq was large enough in magnitude for wave amplihca- 
tion, it had the wrong sign, as shown in Fig. [^. With 
the wrong sign of fcg, the a particle would actually have 
to have higher concentration on the plasma periphery 
than the plasma center for amplihcation to occur. Thus, 
even when the most favorable possible conditions are im¬ 
posed in terms of the magnitude of the a-particle gra¬ 
dient, namely |Aa| < |Ao|, the wave is strictly damped 
rather than amplified on the a particles. To accentuate 
this point, we show in Fig. how even in the case of an 
unrealistically and unphysically high peak a particle den¬ 
sity, the wave is damped on the a particles in the case of 
high-held launch. Thus, it is the persistence of negative 
ke in launching from the low-held side that necessitates 
high-held-side launch for LH-mediated a channeling. 

Note also that, in the case of amplihcation, there is 
great sensitivity to the number of a particles, since, in 
the linear limit, amplihcation is exponential-like, with the 
amplihcation rate is directly proportional to the number 
of a particles participating. To illustrate this sensitivity. 



FIG. 5. The evolution of kg (blue), \kr\ (red), and k± 
(green) along the poloidal trajectory length s for low-field- 
side launch in an ARG-like equilibrium. Launch parameters 
are ny = —2.15, / = 4.1 GHz. kg becomes negative rather 
than positive. 


consider, for example, a factor-of-three decrease in the a 
particle density for the same ray as in Fig. This re¬ 
sults in effectively extinguishing the amplihcation effect, 
as shown in Fig. On the other hand, a factor-of-three 
increase in the a particle density for this same ray results 
in a factor of 1.6 x 10® increase in wave power, as shown 
in Fig. [^. Of course, for any initially substantive power 
level, the amplihcation exhibited in Fig.[^ is likely to vi¬ 
olate the assumption of linearity in the rf power, since at 
high intensities, both the electron and a-particle distri¬ 
bution functions would respond to the waves. However, 
this example indicates how in principle large effects could 
take place, which for accuracy would necessitate consid¬ 
eration of nonlinear or quasilinear effects. 


IV. SUMMARY AND DISCUSSION 

What has been shown through the simulations pre¬ 
sented here is that LH waves launched from the tokamak 
high-held side, which are already predicted to be advan¬ 
tageous from an engineering standpoint,^^ can experi¬ 
ence dramatic amplihcation as they penetrate near the 
plasma center. The convectively amplihed wave eventu¬ 
ally damps on electrons, but in the process it takes en- 
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FIG. 6. Low-field-side launch simulation results near the 
plasma center for an ARG-like equilibrium with ny = —2.15, 
/ = 4.1 GHz. For illustrative purposes only and neglected for 
the purposes of calculating the ray trajectory, the a particle 
distribution has the unrealistically large peak density of 2.18x 
cm“®, and a radial decay length of 1mm. (a) Aa (black) 
and Ao (cyan) vs. s, showing |Aq,| < Ao for the entire wave- 
particle interaction region, (b) Wave power vs. s, showing 
damping of the wave. Color available online. 
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FIG. 7. Effect of factor-of-three decrease (a) and increase 
(b) in a particle concentration on wave amplification for an 
ARC-like equilibrium for ny = —1.65, / = 3.7 GHz, and 
9 — n (inside launch). Note change in scale between (a) and 
(b). Golor available online. 


ergy from the a particles leaving it with more energy 
with which to drive current. Effectively, the current- 
drive efficiency is increased, possibly dramatically, if the 
a-particle gradients are large enough. Thus, the effect 
could prove critical to achieving efficient current drive at 
the plasma core. 

The effect occurs as follows: For high-held side launch, 
the poloidal wavenumber has the proper sign for a chan¬ 
neling. Then, as the LH wave encounters its deepest 
penetration, it circles the plasma core, meaning that its 
radial group velocity necessarily vanishes, leaving the 
full perpendicular wavenumber in the poloidal direction. 
Since for the LH wave, the perpendicular wavenumber 
is far greater than its parallel wavenumber, the poloidal 
wavenumber thus grows large enough to capitalize on ra¬ 
dial gradients in the a-particle distribution. At the same 
time as the radial group velocity vanishes, the poloidal 
group velocity is minimized, so that the ray spends more 
time just at the fortuitous location where the channeling 
effect is thus optimized. The confluence of these circum¬ 
stances suggests that efficient channeling can occur for 
inside launch, provided that the ray encounters a suffi¬ 
ciently steep radial distribution of energetic a particles 
at its point of closest approach to the plasma core. In 
practice, the choice of launch parameters might there¬ 
fore be dictated by the location of the steepest a-particle 
gradient. In the example given, for a reasonable total 
number of a particles, and a steep gradient, the wave 


power damped on electrons could be increased by a fac¬ 
tor of 2.53, with most of that damping occurring near 
the plasma core. 

However, there are several caveats worth noting. 

First, our model assumed that the RF waves were not 
strong enough to significantly modify the a-particle dis¬ 
tribution. While valid in the limit of weak RF, this model 
is insufficient in the main cases of interest, namely when 
the diffusion of a particles by the wave occurs on a time 
scale faster than the slowing down of the a particles. For 
if this were not so, then necessarily most of the a-particle 
energy would have been lost to collisions before the waves 
could capture that energy. To achieve significant savings 
in a reactor, such as to support in full the plasma current 
or to achieve the hot ion mode (where the ion tempera¬ 
ture exceeds the electron temperature), the lower hybrid 
waves would be in a regime in which the linear approx¬ 
imation fails. Similarly, it would not be possible to ap¬ 
ply directly these results to the cyclic operation regime 
of current drive ,1^ where synergies were recognized in 
achieving current drive together with a channeling, again 
because the cyclic operation regime would be far from the 
linear regime where these results are valid. Note also that 
the regime of intense LH waves would also necessitate a 
quasilinear treatment of the electrons as well. 

Second, the wave amplification relies on a favorable 
population inversion, which requires for lower hybrid 
waves that the a particle gradients be on the order of 
centimeters. For preliminary reactor experiments, sus¬ 
tained fusion conditions are only marginally met in the 
plasma center, and not at all nearby, which might indeed 
result in a peaked enough energetic a-particle distribu¬ 
tion. However, in an advanced reactor, steep gradients 
are more difficult, and the a-particle distribution may be 
flattened along the diffusion paths when the waves are 
intense enough that a significant fraction of the alpha- 
particle energy is extracted.^ In that case, the a parti¬ 
cles would not slow down significantly through collisions 
before being diffused to lower energy by the LH waves. 
To maintain the sharp spatial gradients, therefore, a sec¬ 
ond wave might be used in addition to the LH wave.^^ 
An example of a second wave that might be useful is 
one that pushes the a particles further in space for the 
same amo unt o f energy extracted, such as the ion Bern¬ 
stein wave j2dl^ possibly even together with waves of even 
lower frequency such as toroidal Alfven eigenmodes.^^. 
Note that the use here of a second wave is essentially dif¬ 
ferent from other uses of a second wave to optimize the 
current drive efficiency,l^^^^^ since, rather than promoting 
absorption by higher velocity electrons, the second wave 
facilitates the absorption of energy to the LH wave from 
the a particles. Both methods should lead to an effective 
increase in the efficiency. 

In sum, although there remain several caveats in apply¬ 
ing these results to the most useful regimes, namely those 
in which the a-particle amplification leads to very signif¬ 
icant channeling of energy, the linear regime addressed 
here is strongly indicative of the potential inherent in 
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LH waves launched from the high-held side. 
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